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Abstract Blooms of the harmful alga, Karenia brevis on the west Florida continental shelf are thought to 
initiate offshore before manifesting as a nuisance along the coastline. Contributing to such blooms are a 
complex sequence of events occurring within oligotrophic waters, which in any given year may or may not 
be facilitated by the ocean circulation. Once initiation occurs, the delivery from the region of offshore 
origination to the region of coastline manifestation requires an upwelling circulation, whereby K. brevis cells 
are advected shoreward along the bottom. The 2018 K. brevis bloom was particularly intense owing to cells 
from the preceding 2017 bloom being reinforced by a newly formed bloom in 2018, a year when the offshore 
conditions in spring through early summer were again favorable for bloom development. As an event 
response to determine the potential for new cells to be delivered to the shore, a glider was deployed from 24 
August 2018 to 17 September 2018 with a track line designed to map water properties over the hypothesized 
initiation region. The coastal ocean circulation during the deployment interval was generally upwelling 
favorable, but the passage of Tropical Storm Gordon temporarily disrupted this flow, after which K. brevis 
appeared along the Florida Panhandle coast. Strong upwelling then reestablished and K. brevis was 
subsequently observed along Florida's east coast. We describe the glider deployment, the K. brevis 
observations, and we use a numerical circulation model to account for the K. brevis manifestation as 
occurred along Florida's west. Panhandle, and east coasts. 

Plain Language Summary We account for the intensity and location of the 2018 Karenia brevis 
red tide outbreak on the west Florida continental shelf by a combination of water property observations and 
numerical circulation model simulations. These confirm the initiation region being offshore, the 
manifestation region being along the shoreline and the delivery mechanism (from initiation to 
manifestation) being an upwelling favorable coastal ocean circulation. The intensity is attributed to the cells 
remaining in the manifestation region from the prior 2017 bloom being reinforced by cells newly formed 
offshore in 2018. 


1. Introduction 

The west coast of Florida experiences blooms of the harmful alga, Karenia brevis, a dinoflagellate whose tox¬ 
ins cause human respiratory distress, kill fish and other sea life, and disrupt the economies of the coastal 
regions where K. brevis blooms manifest. Not all years are subject to similar concentrations nor durations 
of K brevis red tides. For a canonical red tide year, the expectation is for an onset in late summer to early fall, 
with concentrations diminishing into winter. 

The years 1999, 2000, and 2001 behaved in this canonical manner, but most years behave much more ran¬ 
domly. For instance, 2010 experienced virtually no red tide from the spring season onward, 2011 and 2012 
had substantial K. brevis cell concentrations, whereas 2013 showed only a very short-lived event. From 
2014 onward, K. brevis blooms have occurred in each of these subsequent years. It is now fair to say that 
2018 bloom (or the extended bloom that occurred from September 2017 to January 2019) was the worst of 
the red tide events since 2005, when high cell concentrations lasted from winter through summer until 
Hurricane Katrina ushered the beginning of its demise. 

The epicenter for Florida red tide bloom manifestation (Figure 1) is the shoreline region between the Tampa 
Bay and Charlotte Harbor estuaries, although blooms often show up on the Panhandle coast and in some 
years they even appear on Florida's east coast and as far north as the Carolinas. The 2018 bloom is one for 
which K. brevis appeared within the epicenter region, the Panhandle, and the east coast. 
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Figure 1 . A composite of all observations of Karenia brevis made by the Florida Wildlife Research Institute from 1953 
through 2007. Note the epicenter region from just north of the Tampa Bay to just south of Charlotte Harbor estuaries. 
The ellipse denotes the hypothesized, midshelf initiation region, and the arrows denote the transport pathways between 
initiation and manifestation within the epicenter region and along the Florida Panhandle. 


Following a hypothesis on bloom generation, we recently mounted a glider deployment to sample water 
properties over the region thought to be the formative one for Florida red tide and spanning the latitudinal 
extent from which such red tide cells would be transported from the formative region offshore to the man¬ 
ifestation epicenter region alone the shoreline. The deployment interval (24 August 2018 to 17 September 
2018) was at a time when the West Florida Continental Shelf (WFS) was under a state of upwelling, induced 
by the Gulf of Mexico Loop Current interacting with the shelf slope near the Dry Tortugas. The passage of 
Tropical Storm Gordon temporarily disrupted this circulation about 1 week into the deployment, after which 
the circulation returned to its upwelling favorable state. This deployment also coincided with the spatially 
extensive southwest Florida K. brevis bloom that was first observed within the Charlotte Harbor region in fall 
2017, presenting an opportunity to investigate the potential for bloom intensification over the epicenter 
region at the time when new initiation typically occurs. The deployment also spanned the time frame when 
K. brevis concentrations first appeared along Florida's Panhandle coast and subsequently along Florida's 
east coast. 

The present paper describes the water properties observed during the glider deployment and applies a 
Lagrangian particle trajectory analysis to account for the distribution of K. brevis as observed across the west 
Florida epicenter region, along Florida's Panhandle coast and along Florida's east coast. Section 2 provides 
the appropriate background information, section 3 describes the observations, section 4 discusses the trajec¬ 
tories, and section 5 summarizes and makes recommendations for improving K. brevis forecasting. 


2. Background 

Prior to 1998, it was postulated that blooms of Gymnodinium breve, identified as the organism responsible for 
WFS red tide blooms, were initiated some distance offshore (Steidinger, 1975) before appearing in medium 
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(10 5 -10 6 cells per liter) to high (>10 6 cells per liter) concentrations nearshore in late summer to early fall, 
primarily between the Tampa Bay and Charlotte Harbor estuaries. Interdisciplinary studies coincident with 
anomalous upwelling conditions of 1998 led to a fuller understanding of the biological and physical connec¬ 
tions responsible for such blooms. Two companion papers (Walsh et al., 2003; Weisberg & He, 2003) 
explained some of these. The first showed that anomalous upwelling carried cold, deeper ocean water across 
the shelf break and to the nearshore within the bottom Ekman layer. By examining historical 
temperature/nutrient correspondences, it was further noted that these upwelled, deeper ocean waters ven¬ 
tilated the WFS with new inorganic nutrients. The second showed how the new inorganic nutrients were 
consumed across the WFS, thereby favoring faster-growing phytoplankton (e.g., diatoms) over dinoflagel- 
lates and inhibiting the occurrence of a major K. brevis bloom in that year (we note that 
Gymnodinium breve was renamed Karenia brevis in 2001 in honor of the pioneering work by Dr. 
Karen Steidinger). 

The biology of K. brevis is in itself quite complex. How is it possible for a slow growing dinoflagellate to out- 
compete faster growing diatoms? This was hypothesized (Walsh et al., 2006) to occur via a complex sequence 
of events, with the nitrogen fixing cyanobacterium Trichodesmium (helped through iron fertilization by 
Saharan dust, Lenes et al., 2001) feeding K. brevis offshore, where the lack of dissolved silicate precluded dia¬ 
tom dominance. Once established as a monospecific bloom, K. brevis is capable of using diverse sources to 
satisfy continued nutrient requirements, including the use of toxins to kill fish (e.g., Heil et al., 2014; 
Walsh et al., 2009). Thus, an established K. brevis bloom is difficult to extinguish. If transported to the coast¬ 
line via upwelling, the bloom can then utilize riverine color dissolved organic matter for shading and gain 
additional nutrient support from terrestrial or oceanic sources (e.g., Hu et al., 2006; Vargo et al., 2008), con¬ 
tinuing in a bloom state for months. Additional information on the biology of K. brevis are available in two 
special harmful algal bloom (HAB) issues, one by Walsh and Kirkpatrick (2008) and another by O'Neil and 
Heil (2014), and the utilization of satellite ocean color information for characterizing bloom evolution is dis¬ 
cussed by Soto et al. (2018). 

The K. brevis bloom of 2005 was among the most severe events recorded. The bloom onset was offshore, and 
the pathway to the coast was the bottom Ekman layer as demonstrated through model 
simulation/observation comparisons (Weisberg et al., 2009). These results provided further explanation on 
why K. brevis blooms on the WFS manifest between Tampa Bay and Charlotte Harbor, that is, the primary 
destination following the upwelling transport pathway to the coastline. 

Regardless of the biological complexity, it is now clear that the ocean circulation physics play an important 
role in K. brevis ecology because the circulation is what determines whether or not the water properties are 
conducive for the biological complexity to function effectively. For instance, no red tide was observed on the 
WFS in 2010 owing to the anomalously strong and protracted upwelling that occurred in 2010 (Weisberg, 
Zheng, Liu, et al., 2014). What Weisberg and He (2003) and Walsh et al. (2003) jointly discussed for 1998 
occurred again in 2010, and given improved modeling tools (Zheng & Weisberg, 2012), we were able to make 
the case for why no bloom occurred. With these concepts established, compelling arguments were offered for 
why the K. brevis bloom of 2012 was robust, whereas the one in 2013 was relatively nominal and short lived 
(Weisberg, Zheng, Liu, Corcoran, et al., 2016). These coastal ocean circulation arguments were then 
employed successfully (albeit unpublished) to predict the 2014 and 2015 K. brevis blooms. Seeking a more 
formal method for seasonal prediction Liu et al. (2016) implemented a neural network analysis, combining 
satellite altimetry data with historical K. brevis cell counts for the period 1993-2015. To date, this approach 
accounted for either the occurrence or lack of occurrence of major K. brevis blooms for 20 of the past 25 years. 

Whereas the region from Tampa Bay to Charlotte Harbor appears to be the epicenter for K. brevis HABs on 
the WFS, there are instances when we see blooms along Florida's Panhandle coast and even along Florida's 
east coast. Both the Panhandle and east coast blooms may share their origin with the epicenter region 
blooms. A period of northward advection can transport cells from the northern portion of the initiation zone 
to the Panhandle coast, and a period of protracted upwelling can transport cells around the Florida Keys and 
into the Gulf Stream. Once in the Gulf Stream, these cells can arrive at the east coast within a very short time 
duration (only a few days to Palm Beach and only a week or two to the Carolinas). Examples of Florida 
Panhandle blooms, also with a bottom Ekman layer delivery to the nearshore, are given by Kamykowski 
et al. (2013) and McCulloch et al. (2013), and an example of a Florida east coast bloom, emanating from 
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Figure 2. K. brevis cell counts plotted using log scale from 1 January 2017 through 16 December 2018 (from Florida Fish 
and Wildlife Conservation Commission - Fish and Wildlife Research Institute HAB monitoring database) show the 
persistence of the K. brevis bloom from 2017 through the spring and early summer of 2018 in the epicenter region 
(black circles), and the initiation of bloom conditions in the Panhandle (red) and on the east coast (blue). Only samples 
with >10 3 cells per liter are shown (data source: Florida Fish and Wildlife Conservation Commission - Fish and Wildlife 
Research Institute HAB monitoring database). 

the WFS epicenter region, and traveling around the Florida Keys to as far north as the Carolinas, is given by 
Walsh et al. (2009). 


3. Observations 

With the exception of 2013, major WFS K. brevis blooms (sustained concentrations exceeding 10 6 cells per 
liter) have occurred annually since 2011. The 2017 bloom started in late summer/fall as anticipated, but 
for reasons not understood, it persisted throughout the winter and subsequent spring (Figure 2). With the 
Loop Current penetrating well into the Gulf of Mexico and distant from the Dry Tortugas region during 
the spring and early summer months of 2018, it was surmised (based on Liu et al., 2016) that 2018 would also 
experience a major K. brevis bloom, perhaps rivaling that of 2005 because whatever grew anew for 2018 
would be additive to what was already in place from 2017. 

To explore this possibility, we mounted the deployment of a profiling underwater glider, hereafter simply 
referred to as a glider. Gliders profile the water column from surface to near bottom for weeks at a time, tra¬ 
versing 100 s of km (Rudnick et al., 2004). With a suite of sensors, they sample water column variables and 
transmit these data back to the user several times per day. University of South Florida's Teledyne Webb 
Research Slocum glider was deployed from 24 August 2018 to 17 September 2018, with the intent of mapping 
the region offshore from which it was thought that new cells would be transported to the nearshore epicenter 
region. The deployment started at about the 25-m isobath with a track line planned to zigzag from the initial 
position just to the north of Tampa Bay southward to the latitude of Sarasota, Florida. Figure 3 shows the 
glider track, which deviated somewhat from the original plan owing to the passage of Tropical Storm 
Gordon on 4 September 2018. 

Prior to the glider deployment, the Loop Current shed an eddy in mid-July and retreated back to its more 
direct inflow (through the Yucatan Strait) to outflow (through the Straits of Florida) route, which placed 
it in contact with the Dry Tortugas region “pressure point.” While this contact came too late in the summer 
to suppress a new 2018 K. brevis bloom, its effect by setting the WFS in an upwelling circulation state ensured 
that whatever was blooming along the bottom offshore would eventually appear along the coast. This upwel¬ 
ling effect was observed throughout the month of September 2018, as shown in the Figure 4 composite of 
K. brevis cell concentrations. What was previously observed to the south of Tampa Bay and in particular 
to the south of Venice as holdovers from the 2017 bloom came into full effect over the entire epicenter region. 
Moreover, high cell concentrations also appeared along the Panhandle coast after Tropical Storm Gordon in 
early September and then on Florida's east coast toward the end of September (Figure 3). 

The glider sampled temperature, salinity, chlorophyll, and dissolved oxygen, for which transect plots are 
provided in Figure 5. The low temperatures and higher salinities observed at depth are a reflection of the 
cooler, saltier waters transported across the shelf within the lower portion of the water column by the upwel¬ 
ling circulation. The highest values of chlorophyll are also observed near the bottom along with the lowest 
values for dissolved oxygen. While increased chlorophyll values are not uncommon at depth, here it is noted 
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Figure 3. The glider mission track from 24 August 2018 to 17 September 2018. The loop observed around 4 September is 
attributed to tropical storm Gordon whose resulting northward currents were too strong for the glider to compete against. 

that the highest values do not appear to be contiguous with the deeper ocean; instead these values are 
roughly centered between the 30- and 40-m isobaths, that is, they appear to be of midshelf origin. Where 
the chlorophyll is high, the oxygen tends to be low, indicative of respiration and/or the oxidation of 
decaying organic matter. While these factors alone do not confirm a K. brevis origin, a limited number of 
water samples taken during the deployment confirm their K. brevis nature, and while these limited 
samples are insufficient for a detailed calibration determination, they do suggest that the glider was 
observing concentrations in excess of 10 4 cells per liter. Also of interest are the diurnal variations in 
dissolved oxygen with highest values in daytime and the diurnal variations in chlorophyll indicative of 
vertical migration with cells leaving the bottom during hours of darkness (e.g., Heil, 1986; Hu et al., 2016). 

4. Interpretations 

The observation of near-bottom chlorophyll maxima, later identified through microscopy as associated with 
K. brevis (at least for the limited number of water samples that were obtained) support the hypothesis of the 
K. brevis formative region being at midshelf and displaced northward from the manifestation epicenter 
region located primarily between the Tampa Bay and Charlotte Harbor estuaries. Based on the circulation 
being the delivery mechanism from formation to manifestation, can we account for the composite of 
September 2018 K. brevis observations shown in Figure 3? To answer this question we employ the West 
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Figure 4. K brevis cell count composite for the month of September 2018. 


Florida Coastal Ocean Model (WFCOM), which consists of the unstructured grid, Finite Volume 
Community Ocean Model (FVCOM, e.g., Chen et al., 2003) nested in the Gulf of Mexico Hybrid 
Coordinate Ocean Model (GOM-HYCOM, e.g., Chassignet et al., 2009). FVCOM, as an unstructured grid 
model with flooding and drying, affords very high resolution and accommodates both realistic coastlines 
and shallow water depths. Thus, it allows for inclusion of the estuaries and Florida Bay, along with the inlets 
connecting these with the Gulf of Mexico and across the Florida Keys to the Straits of Florida. GOM- 
HYCOM, as a data assimilative, structured grid model inclusive of tides, facilitates the inclusion of deeper 
ocean forcing for the WFS. By nesting FVCOM in GOM-HYCOM, we downscale from the deep ocean, across 
the continental shelf and into the estuaries. Zheng and Weisberg (2012) provided a WFCOM proof of con¬ 
cept, and subsequent applications include gag grouper recruitment (Weisberg, Zheng & Peebles, 2014), 
K. brevis red tide (Weisberg, Zheng, Liu, et al., 2014; Weisberg, Zheng, Liu, Corcoran, et al., 2016) and track¬ 
ing of the Deepwater Horizon oil spill (Weisberg, Zheng, Liu, Murawski, et al., 2016; Weisberg et al., 2017). 

WFCOM presently provides daily, automated nowcasts and forecasts for the coastal ocean circulation from 
west of the Mississippi River to south of the Florida Keys. These also include short-term (4.5-day) 
hindcast/forecasts of K. brevis locations initialized with observations provided by the Florida Fish and 
Wildlife Conservation Commission-Fish and Wildlife Research Institute as a product of our Collaboration 
on the Prediction of Red tide. Here we employ WFCOM in hindcast to determine if we can account for the 
Figure 4 observations. Lagrangian trajectories are calculated using the isopycnic particle trajectory formal¬ 
ism of Weisberg, Zheng, and Liu (2016) that constrains particles to follow surfaces of constant potential den¬ 
sity. The isopycnic Lagrangian calculations are performed off-line using hourly averages saved from the 
Eulerian WFCOM run. Thus, while mixing may alter the density field throughout the simulation, for each 
of the Lagrangian (hourly) time steps, the kinematic constraint of no flow across isopycnals is adhered to. 

Results for particles initialized at the near-bottom level (ct = —0.9) at locations and times coinciding with the 
glider deployment are shown in Figure 6. Recall that the glider deployment goal was to map out the 
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Figure 5. From top to bottom are temperature, salinity, chlorophyll, and dissolved oxygen, as observed along the glider 
track from 24 August 2018 to 17 September 2018. Note that (1) the salinity was corrected for the conductivity sensor 
time constant using the procedure provided by Liu et al. (2015) and (2) the line on 10 September coincides with a cleaning 
of the chlorophyll optics. 


hypothesized K. brevis formative region from which K. brevis cells might be advected to the manifestation 
epicenter region along the shoreline. Based on the circulation over the diagnosed time interval, the 
hypothesis on transport from initiation to manifestation is demonstrated to be correct. 

As a second hypothesis test, we may ask whether or not cells from the hypothesized formative region could 
have reached the Florida Panhandle coast. Whereas the protracted upwelling conditions continued along 
the bottom in the model simulation, the surface currents did reverse for several days beginning prior to 
the passage of Tropical Storm Gordon, and with vertical migration, it is plausible that cell concentrations 
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Figure 6. Isopycnic trajectories for particles initialized along the glider track from 24 August 2018 to 17 September 2018 
and run through 1 October 2018, a total of 2 to 5 weeks. The initial vertical location for all of these particles was the near¬ 
bottom (a = —0.9) level, and the daily color coding provides the actual depth en route. 


of lower magnitude may have occurred throughout the water column. Assuming this to be correct, Figure 7 
demonstrates that cells migrating to the surface could have been advected to the Panhandle over the period 
associated with Tropical Storm Gordon, assuming that the formative region extended somewhat farther 
north than that covered by the glider deployment. 

As a final hypothesis test, we may ask if the bloom subsequently seen on Florida's east coast could also have 
arrived from the formative region. Two scenarios are possible. The first assumes a direct pathway from the 
formative region to the east coast; the second assumes an indirect path such that particles first arrived at the 
manifestation epicenter region, later to be transported both offshore and around the Florida Keys. Both sce¬ 
narios are plausible. While not shown, particle tracking results similar to Figure 6, but for cells initialized at 
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Figure 7. Isopycnic trajectories for particles initialized at the surface layer on 25 August 2018 (left-hand panel) or 26 
August 2018 (right-hand panel) and tracked through the end of September. Only those particles released across the 
northern portion of the Big Bend region made it to the Panhandle. 
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Figure 8. Isopycnic trajectories for particles initialized at the surface layer and tracked for 1 month. From left to right are initializations beginning on 1 July 2018,1 
August 2018, and 1 September 2018, the latter being the most plausible of these. 
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Figure 9. Modified Moderate Resolution Imaging Spectroradiometer (MODIS) fluorescence line height (or algal bloom 
index, Hu & Feng, 2016) images collected on 28 August 2018 (a) and 18 September 2018 (b) showing K. brevis extending 
over the bloom manifestation region. Overlaid are K. brevis concentrations (cells per liter) determined by FWRI using 
water samples collected during the same week of the MODIS image collection, which confirmed that ABI values of >0.02 
correspond to >10 3 cells per liter of K. brevis. Note the increase in offshore extent on 18 September relative to 28 August 
owing to strong easterly winds. The inset in (b) is the same image on 18 September but with coverage further to the south 
showing the surface transport of K. brevis cells to the Florida Keys, consistent with the model simulation results in 
Figure 8c. 

midwater column, versus the near-bottom, exhibited advection to the Florida Keys. Figure 8 shows results 
for the indirect scenario. In particular, there was a period of time in early September when the winds 
were persistently easterly, thereby driving the nearshore cells offshore at the surface (Figure 9) only to be 
replenished by new cells arriving along the bottom. Once advected offshore, these near surface particles 
(cells) were entrained in the (Loop Current-induced) southward-directed upwelling flow and advected 
around the Florida Keys. 


5. Summary and Conclusions 

The 2018 K. brevis HAB on the west coast of Florida was perhaps the worst red tide occurrence since 2005. 
We explain this as the confluence of cells remaining from the former 2017 bloom being reinforced to by a 
newly formed bloom in 2018. Such blooms are thought to initiate offshore at midshelf when the water con¬ 
ditions there are nutrient deplete (oligotrophic), allowing K. brevis, through a complex sequence of events, to 
out-complete faster growing phytoplankton. The water conditions in which this may (or may not) occur are 
determined by the ocean circulation. If conducive for bloom development to occur, then there will likely be a 
major red tide event; if not conducive, then there likely will not be a major red tide event in any given year. 
The circulation was conducive in both 2017 and 2018, and hence, we experienced major red tides in each of 
these years continuing through January 2019. Hindcast analyses show that such hypothesis accounts for 20 
out of 25 years for which simultaneous K. brevis and ocean circulation data exist. These conditions are deter¬ 
mined by upwelling, a requirement for transporting new cells from their initiation region offshore to their 
manifestation region along the shoreline. The process is subtle, however, because too much upwelling in 
spring to early summer months, when new K. brevis blooms tend to form, will suppress such formation by 
the advection of new inorganic nutrients onto the shelf from the deeper ocean, thereby favoring faster grow¬ 
ing, benign phytoplankton over the harmful K. brevis. 

As a test for offshore initiation and confluence between blooms from two different years, a glider was 
deployed to map water properties offshore within the hypothesized formative region. As expected, 
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relatively high chlorophyll and low oxygen levels were observed along the bottom at midshelf, along with an 
upwelling favorable circulation. Shortly thereafter, elevated levels of K. brevis were detected along the 
Pinellas County, FL, shoreline (just north of Tampa Bay). A Lagrangian trajectory analysis, using a numer¬ 
ical coastal ocean circulation model, demonstrated the connection (via upwelling) between the offshore gen¬ 
eration region and the coastline manifestation region. Moreover, such analyses also accounted for the 
subsequent appearance of K. brevis red tide, both along Florida's Panhandle and east coasts. Advection to 
the Panhandle may be attributed to Tropical Storm Gordon, and advection to the east coast may be attribu¬ 
ted to continued upwelling and subsequent transport around the Florida Keys and into the Gulf Stream. 
Once in the Gulf Stream, it is only a matter of days before cells can arrive at Palm Beach, Florida, where 
the 2018 bloom was first observed on the east coast, and points farther north. Thus, we can explain why 
2018 was one of those rare years in which red tides of a common origin were simultaneously observed in 
all three of these regions: the west Florida epicenter from Tampa Bay to Charlotte Harbor, Florida's 
Panhandle coast, and Florida's east coast. 

Further advances in our ability to forecast K. brevis blooms for Florida (and elsewhere) will require contin¬ 
ued observations, particularly of water properties across the water column within the formative region. 
Repeated glider sections are an effective way of doing this, along with nutrient sensors mounted near bottom 
on strategically located moorings. From a biological perspective, missing are testable hypotheses on eventual 
bloom demise. Some of this may also be physical. If the convergence of cells both by a combination of advec¬ 
tion and growth/mortality is positive, then a bloom will continue; if negative, then a bloom will abate. Only 
though well-conceived, interdisciplinary programs will answers to these and other questions of coastal ocean 
ecology be obtained. 
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